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Abstract 

The problem of center-of-mass (CM) contaminations in ab initio nuclear structure calculations using configuration interaction (CI) 
and coupled-cluster (CC) approaches is analyzed. A rigorous and quantitative scheme for diagnosing the CM contamination of 
intrinsic observables is proposed and applied to ground-state calculations for "^He and '^O. The CI and CC calculations for '^O 
based on model spaces defined via a truncation of the single-particle basis lead to sizable CM contaminations, while the importance- 
jtruncated no-core shell model based on the Nimj^hQ. space is virtually free of CM contaminations. 
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,1. Introduction 

■ In the realm of quantum many-body systems the atomic nu- 
cleus poses unique challenges that do not appear in the theoret- 
. ical description of other systems. One of them results from the 
fact that the nucleus is a finite self-bound system. Unlike, e.g., 
the many-electron systems in atomic, molecular, and condensed 
matter physics, and chemistry, in which the binding originates 
from an external potential created by much heavier atomic nu- 
clei that define a reference frame, the nucleus is bound solely 
by the interactions between the constituent nucleons without 
any external confinement. 

In the exact theory, the many-body state I*!*) of the nucleus 
Ifactorizes into the intrinsic state |i^int) and a state |i^cm) describ- 
ing the dynamics of the center-of-mass (CM), 



X 



m = liAin.) ® l-Acm) 



(1) 



Dnly this complete decoupling guarantees that the intrinsic state 
of the nucleus is translationally invariant and, as result, that all 
intrinsic observables are free of spuriousities induced by the 
CM state of the system. The dynamics of li/'im) is governed 
solely by the intrinsic Hamiltonian 



(T- 



Tcm) + V ^Ti„, + V, (2) 
^ is the intrinsic kinetic energy 



where Jim = jijjx i:,<>(P/ - P;) 
and V is the nuclear interaction. 

Generally, the wave function factorization given by Eq. ([1]) 
no longer holds in approximate calculations. One way to en- 
force it is through the use of Jacobi coordinates, as in, e.g.. 
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the translationally invariant formulation of the no-core shell 
model (NCSM) lll|2[], the hyperspherical harmonics approach 
Istl, and, implicitly, also the Green's function Monte Carlo 
method 101 ■ Due to the computational complexity, these ap- 
proaches are limited to small few-body systems. For this rea- 
son, most approaches aimed at heavier nuclei adopt some finite- 
dimensional model space spanned by Slater determinants con- 
structed with a single-particle basis set, such as the harmonic 
oscillator (HO) basis. Generally, the individual Slater determi- 
nants are not translationally invariant and most model spaces 
do not restore this invariance automatically, i.e., the solution of 
the many-body problem for Hi^t leads to states with a coupling 
between intrinsic and CM motions which is induced by and de- 
pends on the structure of the model space. The one exception 
is the A^niax^f^ space of NCSM 11112]. For a Slater-determinant 
basis of HO single-particle states truncated with respect to the 
HO excitation energy NhQ. there exists a unitary mapping onto 
a Jacobi-coordinate basis. Thus, the exact factorization of Eq. 
([T]) is possible in the Nma^hQ. space. For other single-particle 
bases or model-space truncations, the factorization is lost and 
a coupling of intrinsic and CM states emerges. In this commu- 
nication, we provide a quantitative analysis of the CM contam- 
ination in the truncated ab initio configuration interaction (CI) 
and coupled-cluster (CC) calculations for "^He and "'O. 

2. Center-of-mass diagnostics 

Because a direct analysis of the CI and CC wave functions 
with respect to the factorization given by Eq. ([TJ is not feasi- 
ble, we need a practical tool to assess the degree of unphysical 
coupling between intrinsic and CM motions resulting from the 
truncations used in the CI and CC calculations. A stringent yet 
simple probe can be proposed based on replacing Hint entering 
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the CI and CC calculations by 



Hp - Hu 



(3) 



which includes a Hamilton-type operator H^m acting exclu- 
sively on the CM part of the many-body state of interest, with 
parameter /3 controlling its strength. In principle, any operator 
depending on the CM position Xcm and momentum Pcm, with a 
spectrum bounded from below, could be used for this purpose, 
but we further require that the exact ground state of //cm can be 
represented in the model space. Thus, in a ground-state calcu- 
lation with Hjs for non-zero /3 the CM is in its exact ground state 
if the corresponding many-body wave function |^^) factorizes. 

For calculations based on Slater determinants of HO single- 
particle states, an operator which meets these requirements is 



Ha 
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-X 



-HQ. . 



(4) 
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as was first adopted by Palumbo fst) as well as Gloeckner and 
Lawson ^ in the context of the CM problem in the valence- 
space shell model 101 ■ In an Nmnxfi^ space, the lowest A^max + 1 
eigenstates (including both parities) of //cm, Eq. (|4]i, are repro- 
duced exactly. Hence, as long as the OhQ space is a subspace of 
the model space used in the CI or CC calculation, the solution 
of the //cm eigenvalue problem leads to the exact CM ground- 
state energy, which is zero by construction. For a closed-shell 
nucleus, the use of the OpQh determinant as a reference state 
guarantees this. 

By solving the Schrodinger equation for Hp, Eq. using 
different f3 values we can determine to what extent the intrin- 
sic and CM components of the many-body wave function of 
interest are coupled. If a given many-body method leads to a 
factorized state as in Eq. ([T]i, the intrinsic component |i/fint) 
and the intrinsic observables become independent of /3. In this 
case, the /3//cm part of Hp can only affect the CM component 
of I*!*), which has no effect on intrinsic properties. However, 
if the many-body scheme used to solve the Schrodinger equa- 
tion for I*!*) cannot factorize intrinsic and CM motions, the in- 
trinsic observables, most notably the expectation value of Hint, 
will acquire an unphysical dependence on /3. To quantify the 
strength of this unphysical coupling and its impact on the in- 
trinsic energy, we need to monitor the expectation value (//mt);3 
computed with the eigenstates \^'p} of Hp. A necessary condi- 
tion for the wave function factorization given by Eq. ([T]i is that 
{Hmt)p is independent of Thus, we adopt the change of the 
intrinsic energy, A{Hm)p = {Hmi)p - (//int>o, when going from 
y6 = to a finite as a primary criterion for assessing the CM 
contamination in the wave function. This criterion provides a 
quantitative measure of the impact of the CM contamination on 
the intrinsic energy, while identifying states I*!*) that factorize 
according to Eq. ([TJ which satisfy A{Hint)p - 0. 

As a secondary criterion we use the expectation value of //cm 
obtained for non-zero values of f3. If a given ground-state cal- 
culation allows for a factorization of the wave function, we can 
expect (//cmV exactly equal to the lowest eigenvalue of 
//cm, i e-, zero in the case of H^m defined by Eq. (|4]i. The non- 
zero {Hcm}p value at a finite /? indicates a coupling of CM and 



intrinsic motions. We must emphasize, however, that the use of 
{Hcm)p alone is treacherous, since it does not provide any quan- 
titative measure of the effect of the coupling of the intrinsic and 
CM motions on intrinsic properties. The expectation value of 
//cm at /? = 0, which is sometimes used to judge the magnitude 
of CM contamination in approximate many-body calculations, 
does not provide any quantitative information, since (//cm)/3=o 
can assume any positive value when the calculated wave func- 
tion 1*1*) is factorizable. 

We use Hp, Eq. (|3]i, solely for probing the spurious CM cou- 
pling at moderate fi values, although its original use as a way 
to eliminate the CM contamination was tied to large values of 
P (~ 10^) ill. In a CI calculation, the use of Hp with large p 
amounts to a projection of the model space onto the largest em- 
bedded Nmayfi^ space. In comparison to NCSM, this approach 
is inefficient and numerically problematic. 



3. Many-body methods 

We employ the above diagnostics to analyze CM contam- 
inations in ab initio CI and CC calculations, which rely on 
the model spaces spanned by Slater determinants constructed 
from the HO single-particle states that satisfy the condition 
e ^ ^max, where e - 2« + / is the principal quantum num- 
ber, i.e., model spaces that violate translational invariance from 
the outset. The first method we consider is the importance- 
truncated CI (IT-CI) approach introduced in Refs. isilgi fioll in 
which the CI eigenvalue problem in the above model space is 
further reduced through the use of an importance measure for 
the individual many-body basis states derived from perturba- 
tion theory. We focus on the IT-CI approach with up to 4p4h 
excitations from the initial reference state lO) [YT-Cl(4p4h)], 
which for the closed-shell nuclei examined in this work is the 
determinant obtained by occupying the A lowest single-particle 
HO states. In the limit of vanishing importance threshold, IT- 
CI(4/:'4/2) becomes equivalent to the CI method with all singly, 
doubly, triply, and quadruply excited determinants relative to 
lO) [CISDTQ s Cl{4p4h)]. In this study, we follow the recipe 
described in detail in Refs. [l^l, in which we perform a se- 
quence of YT-Cl{4p4h) calculations with decreasing thresholds 
and extrapolate the results to the vanishing importance thresh- 
old limit, thus obtaining a very good approximation to a com- 
plete CISDTQ calculation. 

The IT-CI approach offers two important advantages for the 
present study. First, it allows for a computationally inexpen- 
sive evaluation of the exact quantum-mechanical expectation 
values (//int)/3 and (//cm);3 which require much less effort than 
the corresponding CISDTQ calculations. Second, we can apply 
the same importance-truncation methodology to other model 
space truncations, such as the A^^ax^fi space of NCSM. In the 
limit of vanishing importance threshold, the resulting sequential 
importance-tuncation IT-NCSM(seq) scheme, which we make 
use of in this study as well, reproduces the model space of the 
full NCSM approach jl^ . 

The second many-bod y appro ach we employ is the single- 
reference CC theory fTll [l2l[ll Q. The specific CC models 
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considered in this study include CCSD (CC singles and dou- 
bles) ilSll and CR-CC(2, 3) | 16ll. both used in our recent nuclear 
structure work H [13, III uJln the CCSD approach, the clus- 
ter operator T defining the CC ground state I*!*) = exp(r)|cl'), 

where in the exact case T - Tx+Ti-^ H and T„ designates 

the npnh component of T, is truncated at the 2p2h clusters T2, 
i.e., T -T\ +T2- The T[ and T2 clusters are determined by solv- 
ing the system of non-linear, energy-independent equations ob- 
tained by inserting the CC wave function, in which T = T\+T2, 
into the Schrodinger equation and projecting on the singly and 
doubly excited determinants relative to lO). The CCSD energy 
is calculated afterwards using the converged T\ and T2 clus- 
ters. Since clusters are important to obtain a quantitative 
description of closed-shell systems, we correct the CCSD re- 
sults for their effect through the suitably defined non-iterative 
corrections to CCSD energies, calculated using the CCSD val- 
ues of T\ and T2, defining the CR-CC(2,3) approach [16]. The 
main advantages of CC methods, as compared to truncated CI 
approaches, are size extensivity and the computationally effi- 
cient description of the higher-order components of the wave 
function beyond a CI model space truncated at the same exci- 
tation level via products of cluster operators. For example, the 
basic CCSD method describes ApAh excitations through, e.g., a 
product of two T2 clusters. The details of the CCSD and CR- 
CC(2,3) methods and their comparison to IT-CI can be found in 
Ref. §]. 

In the CC approaches the evaluation of expectation values 
of operators is more involved than in the CI methods. As ex- 
plained in Ref. we have to rely on a response formulation or 
the equivalent derivatives of CC energies, since the traditional 
expectation value expression with the CC wave function pro- 
duces a non-terminating power series in cluster amplitudes that 
does not lead to practical computational schemes. Moreover, 
non-iterative CC methods, such as CR-CC(2,3), rely on cor- 
rections to the energy only and do not have the corresponding 
wave function. Thus, to determine the CC analogs of (//int);3 
and {Hcm}j3, we proceed as follows. First, we evaluate the 
CC energy E/^' for the Hamiltonian Hj3> = Hi„( + /3'Hcm with 
f}' - p - Ay8, /?, and ^ + A/3 around the nominal value of /3 we 
are interested in, where the step equals 0.01. Then, using a 
centered finite-difference form to approximate the first deriva- 
tive, we compute {Hcm)/3 as {dE0> /df3')j3'=0 and from that obtain 
{Hiat)p - Ep - yS(//cm)/3- We have checked the stability of the 
finite-diff'erence calculations of (//int)/s and (//cm)/j by including 
up to five p' points around /? and using a few different A/S values. 
A cross-check of this prescription was also performed for the 
IT-CI approach, where a comparison with the exact expectation 
values was possible. The experience with performing CC calcu- 
lations in quantum chemistry is that the difference between the 
conventional quantum-mechanical expectation values and the 
corresponding energy derivatives, as described above, are very 
small, since the approximate CC methods, such as those used 
in this work, provide results close to full CI and full CI, be- 
ing an exact diagonalization, satisfies the Hellmann-Feynman 
theorem. 

All calculations were performed with the Vucom interaction 
introduced in Refs. and employed in Refs. iilloifioll. 
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Figure 1: Dependence of the intrinsic energy change &(Hmt)/3 (a) and of the 
expectation value <//cm>/3 (b) on /3 in IT-NCSM(seq) (•) and lT-CI(4p4/i) (♦) 
calculations for the ground state of '^O. In IT-Cl(4p4/!), fmax = 5 and t}£l = 30 
MeV. In IT-NCSM(seq), = 8 and ?jn = 30 MeV. 



It is a pure two-body interaction derived from the Argonne 
V18 potential 122] via a unitary transformation to account 
for short-range central and tensor correlations |23]. We con- 
sider a few representative oscillator frequencies from the range 
HQ. = 22 - 38 MeV, which includes the minima of the ground- 
state energy of '^O obtained with the many-body methods and 
single-particle basis sets used in this work (see Ref. |0]). 



4. Illustration of the center-of-mass diagnostics 

As a first demonstration of the impact of the model space 
truncation on the coupling of intrinsic and CM motions, we an- 
alyze the IT-NCSM(seq) and n-Cl(4p4h) results for the ground 
state of "'O. The IT-NCSM(seq) approach is based on the 
Nma^hQ. space of full NCSM, which allows for the exact wave 
function factorization given by Eq. ([T]i, but the importance trun- 
cation formally breaks this factorization property and so we 
have to check the magnitude of the resulting CM contamina- 
tion. In the case of IT-CI(4/94/!), the underlying CISDTQ model 
space violates translational invariance even before the impor- 
tance truncation is invoked and so we expect the unphysical 
coupling between intrinsic and CM degrees of freedom to be 
more severe than in the IT-NCSM(seq) case. 

We study the dependence of the expectation values of the 
intrinsic Hamiltonian, {Hnn}/3, and of the CM Hamiltonian, 
{Hcm)i3, on the parameter fi entering Hp, Eq. (O. Figure [T] 
shows this dependence for an IT-NCSM(seq) calculation with 
A^max - 8 and for an \T-C\{ApAh) calculation with Cmax = 5, 
both at fixed hQ. = 30 MeV. For IT-NCSM(seq), the change of 
the intrinsic energy A{Hint)p - {Hmifp - {Hmt)a is always below 
150 keV. Considering the error bars resulting from the thresh- 
old extrapolation, this is consistent with the exact NCSM value 
of zero. The lack of dependence of {Hnn}p on /3 shows that the 
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IT-NCSM(seq) eigenstates satisfy the wave function factoriza- 
tion given by Eq. ([1]) to a very good approximation. In contrast, 
the YT-Cl(4iAh) calculations show a sizable change in (i/int)/s, 
of up to 3 MeV when /3 increases to 30, which is a signature of 
the significant coupling of the intrinsic and CM motions in the 
IT-Cl(4p4h) wave function. 

The expectation values of the CM Hamiltonian, (HcmV, ex- 
hibit complementary patterns. In the case of the IT-NCSM(seq) 
calculations, they are at the level of 70 keV for fi - Q and drop 
to below 20 keV already for/3 = 1. In contrast, the YV-Cl{ApAh) 
calculations result in the (//cmV values of about 1.5 MeV for 
j6 = which are suppressed to 0.5 MeV and 0.2 MeV for P - \ 
and 10, respectively. One might think that an expectation value 
of Hem of 200 keV is sufficiently small to indicate a CM de- 
couphng. However, despite the smallness of (//cm)/? the intrin- 
sic energy continues to change, by ~ 1 MeV when going from 
P - 10 to y6 = 20. This shows that a value of {Hcm}/3 of a few 
hundred keV is not sufficient to claim that intrinsic and CM 
motions decouple. 

The general /^-dependence of the IT-CI(4p4h) results shown 
in Fig. 1 also holds for the CC calculations. Keeping this sys- 
tematics in mind, we simplify further discussion and consider 
(//int)/3 and (//cm)/3 for fi - Q and 10 only, investigating their 
behavior as functions of mass number, model-space size, and 
oscillator frequency. 



Table 1: Center-of-mass diagnostics for tlie ground state of He using tlie 
VucoM interaction. All energies are in units of MeV. 
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5. Case Study: ''He 

As a first case, we study the systematics of the CM contam- 
ination in the IT-CI(4/74/!ll CCSD, and CR-CC(2,3) calcula- 
tions for the ground state of "'He. Table [T] summarizes the ex- 
pectation values of //im and H^m obtained for /3 - and 10 
for different e^ax-truncated model spaces and different oscilla- 
tor frequencies HQ.. Instead of the intrinsic energy at /? = 10, 
we list the energy change A(//int)/3=io- 

Let us first consider the expectation values of for [} -Q. 
In the case of IT-CI(4p4/!), (//cm)/3=o ranges from 0.9 to 3.6 
MeV without a clear trend with respect to model-space size. 
For CCSD and CR-CC(2,3), the (//cm)/3=o values range from 
2.5 to 9 MeV and 2.3 to 7.6 MeV, respectively. These are large 
values when compared to the intrinsic energies. If one used 
{Hcm)p=o as the sole criterion, as has been done in the CC con- 
text (cf., e.g., I24II ). one would have to conclude that the extent 



of CM contamination is dramatic. However, when considering 
the results for yS = 10 a different picture emerges. For Cmax = 6, 
the change in the intrinsic energy A(//int);3=io is on the order of 
200 keV for IT-CI(4p4/i), 900 keV for CCSD, and 800 keV for 
CR-CC(2,3), implying a sizable CM contamination but not as 
severe as indicated by (i/cm)^=o- Furthermore, A(//int)/j=io de- 
creases with increasing Cmax- This behavior is confirmed by the 
values of (i/cm)/3=io, which are on the order of a few tens of keV 
and decrease with increasing Cmax as well. Both quantities indi- 
cate that CM contaminations are larger for Ml - 38 MeV than 



For IT-CI(4p4/!) the uncertainties due to the threshold extrapolation are 
below 0.05 MeV for {//i,,,) and below 0.01 MeV for 



for T£l - 30 MeV, but they are not as large as one might expect 
based on the fact that "'He is a light nucleus. 

The above observations can be explained in the following 
manner The exact solution of the Schrodinger equation, which 
leads to perfect decoupling of intrinsic and CM degrees of free- 
dom, is approached when the CI or CC calculation allows for 
all possible npnh excitations and when e^s^x approaches 00. "'He 
consists of only four particles, so IT-CI(4/94/!) becomes an ex- 
act theory when the importance truncation threshold vanishes 
and Cniax — > The relatively small values of A{Hint}fi=io and 
{Hcm)j3=\o and their systematic reduction with increasing e^nx 
are, therefore, expected and this is what we observe in our IT- 
Cl(4p4h) calculations in which we extrapolate the results to the 
CISDTQ limit. When an issue of size extensivity is of no im- 
portance, as is the case for a four-particle "'He problem, the CR- 
CR(2,3) approach represents a very good approximation to IT- 
Cl{4p4h) or CISDTQ 1^ and thus the systematic reduction of 
A(//int)/3=io and (//cm)/3=io with increasing e^ax and their rela- 
tively small values, though not on the same level as in the IT- 
Cl{4p4h) case, can be expected and our calculations confirm 
this. 

We can draw three conclusions from the ''He results. First, 
the magnitude of (//cm)/3=o is not a viable criterion for assessing 
the extent of the CM decoupling; in some cases, the {Hcm)/i=o 
values do not even reproduce the qualitative trends. Second, 
the change of the intrinsic energy A(//int)/3=io when going from 
/3 = to 10 provides a direct and robust measure of the degree 
of CM contamination and its effect on the intrinsic energy, with 
{Hcm)i3=\Q offering additional insights. Third, the A(//int)/3 and 
{Hcm)i3 values for nonzero yS show a systematic reduction of the 



Table 2: Center-of-mass diagnostics for the ground state of O using the 
VucoM interaction. All energies are in units of MeV. 
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CM contamination as Cmax °°- This makes sense since in 
the Cmax CO Umit IT-CI(4/?4/z) and CR-CC(2,3) essentially 
recover the full Hilbert space of the four-particle system and 
thus effectively approach the exact and factorized solution of 
the Schrodinger equation. 

6. Case Study: i*0 

In order to examine the CM contamination in the CI and CC 
calculations for heavier nuclei, we analyze '^O. A common 
presumption is that CM contaminations are less severe as A in- 
creases because of a general 1 /A-scaling contained in the CM 
Hamiltonian, Eq. (HI. 

The IT-CI(4/?4/!ll, CCSD, and CR-CC(2,3) results are sum- 



For IT-CI(4p4/i) the extrapolation uncertainties are 0.2-0.4 MeV for {//jm) 



marized in Table |2l We first focus on hD. = 38 MeV, which 
yields the minimum energy in the CR-CC(2,3) calculation for 
^max = 7 Igt], the largest model space we are able to handle in 
the CC case. The expectation values of Hem for /? = are very 
large, reaching about 3 MeV forIT-CI(4/:)4/2) and about 15 MeV 
for CCSD and CR-CC(2,3), with no clear trend as functions of 

^max- 

For a quantitative assessment we again resort to calculations 
using Hp with /3 - 10. When ?iQ - 38 MeV, the typical 
MHm}j3=io values in the larger model spaces are about 2.5 MeV 
for n-Cl(4p4h), 3 MeV for CCSD, and 2 MeV for CR-CC(2,3), 
i.e., all schemes suffer from a sizable CM contamination of the 
intrinsic energies. The energy changes A(//int)/3=io are rather 
stable and do not show a clear trend for the different model 
spaces considered here, except for CR-CC(2,3) which seems to 
show a decrease with increasing Cmax- The expectation values 
of //cm for y8 = 10 and HQ. = 38 MeV are in the 250 - 500 
keV range. For IT-Cl{4p4h) and CCSD they decrease with in- 
creasing Cmax, whereas for CR-CC(2,3) they slightly increase. 
Due to this behavior and our computational limitations we are 
unable to draw definite conclusions on a possible reduction of 
the CM contamination when going to model spaces with very 
large emax- Note, however, that unlike in the '^He case, neither 
n-Cli4p4h) nor CCSD or CR-CC(2,3) recover the full many- 
body Hilbert space for '^O in the e^ax — > °° limit, because 
none of these approaches describes all of the relevant npnh ex- 
citations. Therefore, even when e^ax the factorization of 
the resulting ground states is not possible. 

Concerning the dependence on HQ., we observe a systematic 
increase of A(//int)^=io and (//cm),8=io with increasing HQ. for all 
CI and CC methods and all model-space sizes used in this work. 
Although (//cm)/3=o follows this trend as well, a closer inspec- 
tion of the behavior of (//cm)/3=o in comparison to A(//int)/3=io 
again reveals that the expectation value of //„„ at y6 = cannot 
be used to quantify the degree of CM contamination. For ex- 
ample, in the n-Cl{4p4h) calculations for HQ. - 22 MeV and 
^max = 7 we obtain (//cm)/s=o of less than 400 keV although 
A(//int)/3=io is about 1.9 MeV. We find the virtually identical 
A<//int>/)=io value in the CR-CC(2,3) calculations for HQ. = 30 
MeV and e^ax = 5, but the corresponding (//cm)/3=o of 5.38 
MeV is larger than the aforementioned YT-Cl(4p4h) value of 
{Hcm}/3={) by more than an order of magnitude. 

If we compare the general picture emerging from the IT- 
CI(4;94/!), CCSD, and CR-CC(2,3) calculations for "'O, where 
A(//int)/3=io is typically in the range of 1 - 3 MeV, with the anal- 
ogous results for "^He, where A(//int)/3=io is typically 100 - 900 
keV, then the popular belief that the CM contamination is sup- 
pressed as I /A becomes questionable. There are other mecha- 
nisms affecting the CM contamination which counteract a sim- 
plistic 1/A scaling, the truncation of the many-body model 
space at a fixed excitation level n < A and the accuracy of a 
given many -body method relative to full CI, which accounts for 
all npnh excitations, being some of them. The restriction of 
the model space to up to 4p4h states in IT-CI(4/?4/z) provides 



and below 0.01 MeV for <//cm>- 
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a far less complete approximation to the full Hilbert space in 
the case of '^O than in the '^He case, where YT-Cl{4p4h) is vir- 
tually exact. For CCSD and CR-CC(2,3) the situation is more 
complex, since the truncation of the cluster operator T, which 
enters the exponential ansatz for the wave function, at T2 or T3 
brings a variety of product excitations through the disconnected 
clusters, such as {1/2)TiTt, (l/6)T^, etc., which are beyond the 
4p4h excitation level of CI, but still,"CCSD and CR-CC(2,3) are 
not the exact theories, so the problem remains. For the closed- 
shell nuclei described by pairwise interactions, such as '^O, the 
role of higher-order connecfe^f clusters, such as, e.g., T4, in de- 
scribing the intrinsic energy, is expected to be small 1193, but 
it is not entirely clear what the significance of the higher-order 
cluster components of the wave function neglected in CCSD 
and CR-CC(2,3) on the CM motion and its coupling with the 
intrinsic state of the nucleus is. Therefore, the ability of the 
model space used in the CI and CC calculations to represent 
a factorized eigenstate I*!*) is getting progressively worse with 
increasing A. 

7. Conclusions 

We have formulated simple and rigorous criteria for diag- 
nosing the CM contamination in many-body approaches to the 
nuclear structure problem. They are based on the fact that for 
a factorized many-body state the CM component can be ma- 
nipulated, e.g., by adding a CM contribution /?//cm to the in- 
trinsic Hamiltonian Him, without affecting the intrinsic part. 
The existence of an unphysical coupling between CM and in- 
trinsic motions manifests itself in a dependence of the change 
of the intrinsic energy expectation value A{Hi„i}j3 on since 
A(//int)/3 = for all y6 when the nuclear state factorizes. The 
analysis of A{//int)/3 is our primary diagnostic. A secondary 
measure of the magnitude of the coupling between CM and in- 
trinsic motions is the expectation value of //cm at non-zero /?, 
which should be zero for factorizable states. In contrast to these 
two quantities, the expectation value {Hcm}j3 atp-Q, which is 
sometimes used to examine the degree of CM contamination, 
does not provide a meaningful measure since it can have any 
positive value for a many-body state that factorizes. Further- 
more, we have shown examples where (//cm)/s=o is on the order 
of a few hundred keV when A(//int)/3=io is on the order of a few 
MeV, indicating a sizable CM contamination. 

Using our diagnostics we have demonstrated that the IT- 
NCSM approach, which is based on the Nm-^xfi^ model-space 
cutoff, exhibits a very good approximate factorization with 
MHint)i3 below 150 keV for all /? in spite of the fact that the 
importance truncation used to select the dominant contributions 
to the NCSM wave function formally breaks the rigorous fac- 
torization of the Nmaxfi^ space. Switching from the A^max- to 
the Cmax-truncation of IT-CI(4/94/i) leads to a sizable CM con- 
tamination, with A(//int)/j reaching a few MeV as /3 increases. 
Similarly, the CCSD and CR-CC(2,3) calculations show sizable 
CM contaminations, with A{Hi„t}]3=m typically somewhat larger 
than in the IT-CI(4/94/!) case. With increasing emax, the degree 
of CM contamination decreases for "^He, but one does not ob- 
serve the same clear trend for '^O. This difference can be ex- 



plained by the fixed maximum excitation level of the truncated 
CI and CC calculations used in this work. The exact solution of 
the Schrodinger equation, which leads to perfect decoupling of 
intrinsic and CM degrees of freedom, is approached only when 
the CI or CC calculation allows for all possible npnh excitations 



and when e^ 



00. The fixed excitation level n < A used in 



the CI and CC calculations for prevents this. 

Our findings seem to disagree with a recent claim of Hagen 
et al. 12411 of a factorization of the CC ground states to a very 
good approximation, based on the observation that the expecta- 
tion value (//cm)/j={) for a generalized operator Hem defined for 
an oscillator frequency Q. Q. minimizing (//cmV=o is close 
to zero. For '^O described by the N^LO interaction iS] in a 
Hartree-Fock basis with e^ax ~ 18 they obtain the optimum 
{Hcm}j3=o values in the range of about -0.5 to 1 MeV when 
varying HQ. from 24 to 52 MeV. Aside from concerns about the 
interpretation of negative expectation values for a positive def- 
inite operator and the validity of the generalization of Hem, the 
values of {Hcm}]3=o observed in l24ll are not sufficient to claim 
the approximate factorization of the CC ground states. In our 
view, it is more appropriate to monitor A(//int)^ and {Hcm)j3 at 
finite y6 to assess the degree of coupling of intrinsic and CM 
degrees of freedom in truncated CI and CC calculations. 
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